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SELF-SIMILAR  FLUID  DYNAMIC  LIMITS 
FOR  THE  BROADWELL  SYSTEM 


Abstract.  This  report  discusses  a  new  approach  for  the  resolution  of  the  fluid  dynamic 
limit  for  the  Broadwell  system  of  the  kinetic  theory  of  gases,  appropriate  in  the  case  of 
Riemann,  Maxwellian  data.  Since  the  formal  limiting  system  is  expected  to  have  self¬ 
similar  solutions,  we  are  motivated  to  replace  the  Knudsen  number  e  in  the  Broadwell 
model  so  that  the  resulting  model  admits  self-similar  solutions  in  ^  =  x/t  and  then  let 
s.  go  to  zero.  The  limiting  procedure  is  justified  and  the  resulting  limit  is  a  solution  of 
the  Riemann  problem  for  the  fluid  dynamic  limit  equations.  A  class  of  Riemann  data  for 
which  this  program  can  be  carried  out  is  exhibited.  Furthermore  it  is  shown  that  for  the 
Carleman  model  the  complete  program  can  be  done  successfully  for  arbitrary  Riemann 
data. 


§0.  Introduction 


The  Broadwell  system  of  discrete  kinetic  theory  is  given 


by  the  system  of  partial  differ¬ 


ential  equations 


dh  ^  ^ 

dt  ^  dx 
dt  ^  dx 

dh  .  dh 

dt  ^  dy 

dh  _  ^ 

dt  ^  dy 

dh  ^  dh 

dt  ^""dz 
dt  ^  dz 


(^ihh  +  hh  ~~  2/1/2) » 
<^{hh  +  hh  ~  2/1/2)  ? 
(^ihh  +  hh  ~  2/3/4) » 
<r{hh  +  hh  ~  2/3/4)  ? 
o'ifih  +  hh  ~  2/5/5)  1 
^■(/i/z  +  hh  ~  2/5/5) . 


(0.1) 


The  model  describes  a  gas  of  particles  with  identical  masses  moving  along  three  perpendic¬ 


ular  coordinate  axes  with  the  same  speed  c.  Results  of  a  particular  collision  have  the  same 
probability  and  only  binary  collisions  are  considered.  The  functions  /,  =  fi{x,y,t),  i  = 


1, . . .  ,6  denote  the  densities  of  particles  moving  in  the  six  allowed  directions;  (7/2c  is  the 
cross  section  for  binary  collisions. 


For  flows  which  are  independent  oiy,z  and  for  which  /s  =  /j  =  /s  =  /g  the  above  six 


velocity  Broadwell  model  reduces  to  the  simpler  form 


dt  dx 


\{n-hs^) 


df2_^ 

dt  dx 


(0.2) 


where  for  simplicity  we  have  set  c  =  1  and  ^  £  the  Knudsen  number  or  “mean  free 

path”  of  the  gas.  As  the  “mean  free  path”  is  the  distance  between  successive  collisions,  a 
small  mean  free  path  means  the  gas  becomes  less  rarefled  and  a  “macroscropic”  description 
of  the  gas  based  on  fluid  dynamic  Euler  and/or  Navier-Stokes  equations  should  become 
meaningful. 

The  problem  of  rigorously  passing  to  the  fluid  dynamic  limit  has  a  long  history.  Here  we 
give  a  quick  summary  of  relevant  results.  Additional  references  may  be  found  in  the  book 
of  Cercignani  [4]  for  work  on  the  Boltzmann  equation  and  the  review  paper  of  Platkowski 
and  Illner  [12]  with  regard  to  research  on  discrete  velocity  models  in  the  kinetic  theory  of 
gases. 

First  within  the  realm  of  discrete  velocity  models  the  Ceurleman  model  does  allow  for 
rigorous  passage  to  the  fluid  dynamic  limit.  This  wais  done  in  the  work  of  Kurtz  [11].  But 
as  the  Carleman  model  does  not  conserve  momentum  it  is  perhaps  a  poor  test  case. 

For  the  Broadwell  model  the  basic  result  is  due  to  Broadwell  himself  [1].  One  begins  by 
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rewriting  the  system  (0.2)  as 


^(/i  +  /z  +  4/3)  +  ^(/i  ~  /a)  =  0 , 

^(/i  ~ /a)  +  ^(/i  + /2)  =  0,  (0.3) 

Next  one  mahes  the  ansatz  of  travelling  wave  solutions  fi  =  /i(^),/2  =  /2(^)?  /a  = 
f 3(^)1  ^  where  3  will  be  the  speed  of  the  wave.  Substitution  of  this  ansatz  into 

(0.3)  yields  the  system  of  ordinary  differential  equations 

~Kfi  +  /z  +  4/3)'  +  (/i  —  12)'  =  0 , 


“'S(/i  “  h)'  +  ifi  +  /a)*  =  0, 


(0.4) 


We  pose  downstream  and  upstream  positive,  constant  data  (/i,/2,/3)  — + 
as  6  —*  ±00  which  of  course  is  consistent  with  (0.4)  if  and  only  if  the  data  axe  Maxwellians, 
i.e.  f2  =  Since  5  is  a  constant  (0.4)i,  (0.4)2  may  be  integrated  from  —00  to  9  to 


yield: 


“•s(/i  + /2  +  4/3)  +  (/i  - /z)  =  — ^(/i  4-/2  +4/3)  +  (/i  -72)7 


-■s(/i  -  72)  +  (/i  +  72)  =  -5(71  “  72 )  +  (7i  +  72 ) • 

These  two  equations  determine  f\,  f2  ^  functions  of  f3{9)  emd  s.  Substitution  of  these 
functions  into  (0.4)3  will  yield  an  autonomous  scalar  ordinary  differential  equation  for  fz 
with  precisely  two  equilibrium  points  at  7*  =  (7*7* )^^^-  Since  such  boundary  value 
problems  must  possess  solutions  it  follows  that  a  travelling  wave  solution  exists.  The  value 


of  s  is  found  by  integrating  (0.4)  from  — oo  to  +oo: 

+ ft + + ut  -  ft) 

=  -^(/r +/r +4(/r/2')''0 +(/r  -/D.  (o-^) 

-sUt-ft)  +  Ut  ^ft)  =  -<ft  -ft)  +  (ff  +  ft), 

and  solving  the  system  (0.5)  for  s.  These  are  just  the  Rankine-Hugoniot  jump  conditions. 
Once  the  existence  of  a  travelling  wave  is  established,  we  let  e  —*  0+  and  obtain 


The  limit  function  is  a  weak  solution  of  the  limiting  fluid  dynamic  conservation  laws 

|(/,-/2)  +  ^(/.+/2)=0. 

This  is  because  the  limit  function  is  piecewise  constant  possessing  a  jump  discontinuity 
across  the  shock  x  =  st  and  across  x  =  st  the  limit  function  satisfies  the  jump  condition 
(0.5).  In  fact  the  limit  function  is  a  solution  to  the  Riemann  problem  (0.6)  with  piecewise 


constant  initicd  data 


/,  =/f(i<0),  f,=ft(x>0)-, 


/2=/f(x<0),  /2=/+(l>0). 

Here  we  have  used  the  usual  definition  of  weak  solution:  a  pair  of  bounded  meeisurable 


functions  fi ,  f2  on  (— oo,  oo)  x  [0,  oo)  is  called  o  weak  solution  of  the  Riemann  initizil  value 


problem  (0.6),  (0.7)  provided  that 

j(  J  (/i+/2+4{/i/2)‘^^)^+(/i  -/2)^<ix<i<=0, 
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for  all  C°°  functions  with  compact  support  in  t  >  0,  — oo  <  x  <  oo  and 


^lim^/j(a;,t) 


{ 


X  >  0 
X  <  0 


Introduction  of  the  macroscopic  variables  p  =  fi  +  4(/i/2)^/^  +  /2,  pu  =  /i  —  f 2  shows 
that  the  fluid  limit  equations  (0.6),  (0.7)  can  be  written  in  the  form  of  the  Euler  equations 


dp  ^  ^  „ 

~+-(pu)  =  0, 

§i(pu)  +  £(pg(u))  =  0. 


(0.8) 


with  g(u)  :=  i  12(1  +  —  1  and  Riemann  initial  data 


P  =  P  (x  <  0),  p  =  />+  (x  >  0) ; 


pu  =  p  u~  (x  <  0),  pu  —  p'^u'^  (x  >  0) ; 


(0.9) 


ft, 

P^u^:  =  ft-ft. 

It  should  be  noted  that  the  Boltzmsinn  equation  also  possesses  a  travelling  wave  solution 
for  Maxwellian  data  which  axe  close  (Czdlisch  and  Nicolaenko  [2]).  The  data  in  [2]  must  also 
be  consistent  with  relevant  fluid  dynamic  jump  conditions  which  are  the  Rankine-Hugoniot 
jump  conditions  for  a  shock  wave  in  axi  ideal  fluid. 

In  summary  we  see  for  Riemann  data  satisfying  the  Rankine-Hugoniot  jump  conditions 
associated  with  the  fluid  dynamic  limit  equations  (0.6)  (or  equivalently  (0.8))  passage  to 
the  fluid  dynamic  limit  for  the  Broadwell  model  can  be  achieved  (and  with  a  smallness 
assumption  on  the  variation  of  the  data  for  the  Boltzmann  equation  also). 
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What  can  be  Sciid  regarding  the  fluid  dynamic  limit  for  arbitrary  data  or  for  that  matter 
even  the  more  restricted  case  of  arbitrary  Riemann  data? 

For  the  case  of  smooth  data  Inoue  and  Nishida  [10]  have  shown  that  one  can  pass  to  the 
fluid  dynamic  limit  for  the  Broadwell  system  on  a  sufficiently  small  time  inter\’al  to  yield 
a  smooth  solution  of  the  fluid  dynamic  limit  equations.  They  show  compactness  of  the 
parametrized  sequence  {/f,/!?/! }  satisfying  the  Broadwell  system  in  a  space  that  allows 
passage  to  the  fluid  dynamic  limit.  This  result  does  not  cover  the  case  when  the  solutions 
of  the  limit  equations  have  shocks. 

In  a  work  of  complementary  nature  concerning  again  solutions  without  shocks,  Caflisch 
ajid  Papanicolaou  [3]  show  that  a  given  smooth  solution  to  the  limit  equations  can  be 
approximated  by  a  solution  to  the  Broadwell  system  when  e  is  small.  The  approximation 
program  was  recently  carried  out  at  the  level  of  solutions  with  shocks  by  Xin  [15].  He  shows 
that  given  a  piecewise  smooth  solution  (with  finitely  many  noninteracting  shocks)  of  the 
fluid  dynamic-limit  equations,  there  exist  solutions  to  the  Broadwell  system  which  converge 
asymptotically  to  the  fluid  dynamical  solution  2is  e  — »  0-|-.  Conceptually  the  approximation 
program  presupposes  knowledge  of  a  smooth  solution  in  [3]  or  an  admissible  solution  in  [15] 
to  the  underlying  limit  conserA^tion  laws  (0.6)  and  is  intended  as  a  method  to  solve  the 
Broadwell  system  (0.2)  based  on  solutions  to  (0.6).  By  contraist  the  compactness  method 
attempts  to  construct  solutions  of  (0.6)  eis  limits  of  solutions  of  (0.2). 

In  the  research  discussed  here  we  continue  in  the  spirit  of  the  compactness  issue.  The 
goal  is  to  extend  the  success  of  Broad  well’s  original  travelling  wave  idea  to  more  general 
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Riemann  data,  not  necessarily  consistent  with  the  Rankine-Hugoniot  jump  conditions.  The 

idea  is  based  on  the  following  observation:  Any  system  of  conservation  laws 

dF{U)  dGjU) 
dt  dx  ~  ’ 


U  :  (—00, 00)  X  (0, 00)  — >  IR^  ,  F,G  :  IR^  — *■  IR^,  with  Riemann  data 


X  <  0 
X  >  0 


(0.10) 


must  possess  space-time  dilation  invariance:  For  any  positive  constEint  a  >  0,  the  change 
of  variable  (x,t)  —*  (ax,  at)  preserves  both  the  equations  and  the  initial  data.  Hence 
solutions  of  Riemann  problems  are  expected  to  depend  only  on  the  similarity  variable 
^  =  y,  that  is  U(x,t)  =  U(^). 

For  example,  if  one  was  attempting  to  solve  the  Riemann  problem  for  a  system  of 
conservation  laws  by  the  artificial  viscosity  method 


dF(U)  dG(U)  ^  W 

dt  ^  dx  ^  dx'^ 


(0.11) 


one  might  first  consider  substitution  of  the  ansatz  U(x,t)  =  U(^)  into  (0.11).  But  unfortu¬ 
nately  (0.11)  does  not  possess  space-time  dilationeJ  invariance.  For  this  reason  Dafermos 


[5]  suggested  a  new  type  of  “viscous”  limit  problem 

dF(U)  dG(U)  _.d^U 
dt  di 


(0.12) 


which  does  possess  space-time  dilation  invariance.  Substitution  of  U(x,t)  =  U(^)  into 
(0.12)  yields  the  system  of  ordinary  differential  equations 


-^F(U(0y  +  GiUiOr  =  eU" 


(0.13) 
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ajnd  (0.10)  implies  boundary  conditions 


Ui-oo)  =  U-,  £7(+oo)  =  17+. 


(0.14) 


In  papers  [5],  [6]  Dafermos  and  DiPerna  showed  that  for  N  =  2  a  large  class  of  Rie- 
mann  problems  for  hyperbolic  conservation  laws  may  be  solved  as  limits  of  solutions  of 
(0.13),  (0.14)  as  e  — ♦  0-f-.  The  program  has  been  continued  in  the  work  of  Slemrod  [13], 
Fan  [8],  and  Slemrod  and  Tzavareis  [14]. 

In  the  same  spirit  we  easily  recognize  that  the  Broadwell  system  does  not  possess  space- 
time  dilational  im'ariance.  Hence  we  are  motivated  to  consider  an  artificial  Broadwell 


system 


dt  dx 


(0.15) 


which  does  possess  the  desired  space-time  dilational  invariance.  (Of  course  the  same  obser¬ 
vation  is  true  for  the  Boltzmann  equation  and  any  of  the  standard  discrete  velocity  models 
in  the  kinetic  theory  of  gases.) 

We  make  the  ansatz  fi{x,t)  =  /i(0i  =  /2(0’  =  h{C)i  ^  =  f, 

and  substitute  into  (0.15)  to  obtmn  the  system  of  non-autonomous  ordinary  differential 


equations 


-(f  +  i)/^(?)  =  {/^-/i/2)A. 


(0.16) 


-4/;  =  (/./2-/J)/2£ 
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Henceforth,  we  will  use  the  notation  /  =  (/i,/2,/3  )• 


Since  we  wish 


for  z  ^  0  as  < 


0+  we  impose  boundary  data 


/(-oo)=/-, /(+oo)=/+,i  =  1,2,3  (0.17) 

where  f~  —  (fi  ,  f2  ,  ID  =  (ft  ifi)  Maxwellicin  states:  ft  ft  =  ft^- 

System  (0.16),  (0.17)  is  considerably  harder  to  analyze  than  system  (0.4)  obtained  from 
the  travelling  wave  ansatz.  The  reasons  are  obvious:  (i)  (0.16)  is  non- autonomous  in  the 
similarity  variable  ^  and  (ii)  (0.16)  does  not  possess  any  first  integrals  that  allow  to  reduce 
the  number  of  dependent  variables.  It  does  however  posssess  a  simplification.  Since  /*  are 
equilibria,  we  must  have  /(^)  =  f~  on  (— oo,  —1)  and  /(^)  =  f'^  on  (1,  oo).  The  boundeiry 
conditions  at  ^  =  ±oo  are  replaced  by 

/(-i)  =  r.  /(+i)  =  /'",  (0.18) 

and  (0.16)  need  only  be  considered  on  —  1  <  ^  <  1. 

The  goal  now  is  twofold.  First,  to  construct  solutions  of  the  boundary  value  problem 
(0.16),  (0.18)  for  £  >  0  fixed.  Second,  to  show  that  as  e  — »  0-f-  the  solutions  /'  have 
a  limit  which  is  a  weah  solution  to  the  Riemann  problem  (0.6),  (0.7)  for  the  fluid  limit 
system.  Remarkably  the  second  part  of  this  program  is  easier  than  the  first.  We  exhibit 
here  only  a  class  of  Maxwellian  data  /“,  f'^  for  which  solutions  /®  exist  to  (0.16),  (0.18) 
for  all  £  >  0.  However,  for  any  data  for  which  solutions  /'  exist  for  all  £  >  0,  one 
can  e.xtract  a  convergent  subsequence  /*  — ♦  /  a.e.  in  (—00,00).  The  limiting  function  / 
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is  a  local  M''-Xwellian,  fz  =  (/i/2)^^^  a-e.,  and  a  weak  solution  of  the  Riemann  problem 
(0.6),  (0.7).  The  proof  is  based  on  total  vaxiation  estimates  and  use  of  Helly’s  theorem. 

We  note  that  the  approach  given  here  has  some  resemblance  to  a  recent  paper  of  F. 
Golse  [9].  In  his  paper  Golse  makes  the  self-similar  ansatz  for  the  Broadwell  system 

fj{x,i)  =  FjiO/t  ,  i  =  1,2,3  (0.19) 

where  again  ^  =  f-  Substitution  of  (0.19)  into  (0.2)  yields  the  system  of  ordinary  differ¬ 
ential  equations 

-[(^-i)i^i(or  =  (i"3"--F’ii"2)/£, 

-[«  +  =  (F}  -  F,F2)/e  ,  (0.20) 

-liFziOr  =(F,F2-Fi)/2s, 

which  differs  from  (0.16)  in  the  fact  that  the  left  hand  side  of  (0.16)  has  differentiation 
followed  by  multiplication  while  (0.20)  has  the  reverse.  System  (0.20)  then  has  the  same 
property  as  (0.4)  of  possessing  two  first  integrals.  Golse  exploits  this  property  to  show  the 
existence  of  a  solution  Fj,  Fj,  F3  of  (0.20)  analytic  on  — 1  <  ^  <  1.  The  importance  of  the 
result  is  that  it  displays  explicitly  the  large  time  ^^7)  behavior  of  a  class  of  solutions  to 
the  Broadwell  system  (0.2).  The  solutions  /i,/2,/3  of  course  do  not  possess  space-time 
dilational  invariance  and  do  not  appear  relevant  to  solving  the  Riemann  problem  for  the 
limit  fluid  dynamic  system  (0  6),  (0.7). 

The  paper  is  divided  into  five  section  after  this  one.  Section  1  reformulates  the  non¬ 
linear  boundary  \a.lue  problem  (0.16),  (O.IS)  and  provides  information  on  the  behavior  of 
solutions  to  (0.16),  (0.18).  Section  2  provides  the  main  result  on  the  fluid  dynamic  limit 
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problem:  if  for  fixed  Maxwellian  initial  data  f  ~  f~  (x  <  0),  f  =  /■*"  (x  >  0)  (0.16),  (0.18) 
possesses  a  solution  for  all  values  of  the  Knudsen  number  £  >  0  then  the  sequence  of  so¬ 
lutions  {/'(Ol  (0.16),  (0.18)  possesses  a  subsequence  which  converges  bounded  a.e.  in 

— oo  <  X  <  oo,  /  >  0  to  a  solution  of  the  fluid  dynamic  Riemann  problem  (0.6),  (0.7)  (or 
equivalently  (0.8),  (0.9)).  Section  3  provides  a  class  of  data  for  which  (0.16),  (0.18)  does 
indeed  possess  solutions  for  all  values  of  the  Knudsen  number  £  >  0.  The  proof  of  existence 
for  all  £  >  0  is  based  on  the  contraction  mapping  principle  and  a  continuation  argument. 
Section  4  applies  the  program  originally  described  for  the  Broadwell  system  to  the  simpler 
Carleman  model  of  the  kinetic  theory  of  gases.  For  the  Carleman  model  we  show  that  for 
all  Riemann,  Maxwellian  data  the  fluid  dynamic  limit  can  be  achieved.  Finally  Section 
5  provides  an  appendix  describing  the  behavior  of  non- homogeneous  singular  first  order 
equations;  this  appendix  is  used  in  Sections  3  and  4. 

§1.  The  nonlinear  boundary  value  problem. 

In  this  section  we  derive  various  properties  of  solutions  to  the  nonlinear  boundary  value 
problem  (0.16),  (0.18),  which  for  convenience  we  call  (Ve)- 

Since  the  underlying  differential  equations  are  singular,  it  is  worthwhile  clarifying  what 
is  meant  by  a  solution.  The  weak  form  of  (0.16)  takes  the  form: 


+ 

/  fi{r)dT  = 
^  1 

If 

'  Qifir))dT  , 

(1.1) 

K  +  D/^il; 

+ 

ri7 

/  f2{T)dT  = 

>  * 

-f 

’  Qifir))dT , 

(1.2) 

S  1 

1 

+ 

/  f3{T)dT  = 
Ul 

X 

/  Q{f{r))dT, 

(1.3) 
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for  ^1,  ^2  e  [— 1»  1]>  with  the  notation 


QU)  =  fl-hh  (1.4) 

used  throughout  for  the  collision  operator. 

Definition.  The  triplet  /(^)  =  (/i(Oi  /2(0»  hiO)  defined  for  ^e[— 1,1]  is  a  solution 
of  the  boundary-value  problem  {Vc)  if  fi,f2,  /a  eCf-l,  1]  satisfy  the  integral  relations 
(1.1  -  1.3)  for  6,^2  e[— 1, 1],  and  the  boundary  conditions  fji±l)  =  ,  j  =  1 , 2,  3. 

It  follows  immediately  that  such  a  solution  enjoys  the  properties: 

(i)  /i6CM-1,1),  /26C*(-1,1],  /3ecn(-l,0)u(0,l]), 

(ii)  equations  (0.16)  axe  satisfied  for  — 1  <  ^  <  0  and  0  <  ^  <  1, 

(iii)  the  boundary  conditions  (0.18)  are  satisfied  at  ^  =  — 1  and  ^  =  1. 

Also,  solutions  of  (Vg)  satisfy  the  balance  of  mass,  balance  of  momentum,  and  entropy 
production  equations 


-«/i+/2+4/3)'  +  (/i-/2)'=0,  (1.5) 

-«(/.- A )'  +  (/l  4-/2)' =0,  (1.6) 

~^(/i  +  /z  -}-  4/3  ^nf^y  -|-  (/i  inf  I  —  /2  inf 2) 

=  -i(/l  -  /■  A)(Ml  -  MiA) .  (1.7) 

for  ^e(-l,0)  and  ^€(0,1).  The  last  equation  is  obtained  upon  multiplying  (0.16a)  by 
{inf I  +  1),  (0.16&)  by  {inf2  +  1),  (0.16c)  by  A{inf^  -b  1)  and  adding  the  resulting  identities. 
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The  weak  form  of  (1.5  —  1.7)  reads 


-«/.  +  /2  +  +  (/,  -  +  /  (/,  +  A  +  ih)d(.  =  0 , 

-«/.  -  /2)i|: + (fi + /2)ii: + /'V.  -h)di=o, 

•'ft 

-iUitnh  +  /2«n/2  +  4/3«n/j)||;  +  (/.fn/,  -  /2«n/j)|«;  + 

+  /2M2  +  =  --  f'ul  -  fih)(tnfi  -  inhh)ii 

hi  ^  hi 

for  Cl, ^2  e  [-1,1]. 

The  first  lemma  plays  a  key  role  in  the  analysis,  as  it  will  provide  control  on  the  shapes 
of  solutions  of  {Vg). 

Lemma  1.1.  Let  /  =  (/i,/2,/3)  be  a  continuous  on  [—1, 1]  solution  of  (0.16).  Then 

(i)  QifiO)  not  change  sign  on  the  intervals  (-1,0)  or  (0,1). 

(ii)  Q(/(-i))  =  g(/(o))  =  g(/(+i))  =  0. 

Proof.  From  (0.16)  and  (1.4),  it  follows  that  Q{f)  satisfies  the  differential  equation 

+  ^  +  )  0(/)- 

By  uniqueness  for  ordinary  differential  equations,  if  Q(/(Co))  =  0  at  some  point  Co 
then  Q{f)  vanishes  on  the  whole  interval  (—1,0).  The  same  holds  on  the  interval  (0,1), 
hence  (i). 

In  the  integral  relation  (1.3)  set  C2  =  ^  ^  0  and  Ci  =  0.  Then 


(1.8) 

(1.9) 

(1.10) 


Since  fziO,  QifiO)  are  continuous  we  let  ^  0  and  find  /3(0)  =  /3(0)  +  ^Q(/(0)),  i.e. 

<3(/(0))  =  0.  Similarly,  using  (1.1)  and  (1.2)  we  obtain  Q(f(-1))  =  Q(/(l))  =  0.  | 

For  the  remainder  of  the  section  /  =  (/i,/2,/3)  will  stand  for  a  solution  of  (V^)  defined 
on  [—1, 1]  and  corresponding  to  some  e  >  0  and  positive  Maxwellian  data 

<?(/■)  =  )=  0 ,  /* ,  ,  /,*  >  0 .  (M) 

The  components  of  /  enjoy  the  regularity  :  /i  eC[-l,  l]nC^[-l,  1),  f2  eC[-l,l]nC\-l,l] 
and  fz  e  C[— 1, 1]  fl  C^([— 1, 0)  U  (0, 1]).  Lemma  1.1  together  with  the  form  of  the  equations 
(0.16)  impose  restrictions  on  the  shapes  of  the  functions  allowing  only  for  the  following 
possibilities: 

Ci:  Q{f)  >  0  for  — 1  <  ^  <  0  and  0  <  ^  <  1:  In  this  case  f\  is  increasing  on  (—1, 1), 
is  decreasing  on  (—1, 1),  and  fz  is  decreasing  on  (—1,0)  and  increasing  on  (0, 1). 

C2:  Q{f)  <  0  for  — 1  <  ^  <  0  and  Q{f)  >  0  for  0  <  ^  <  1:  In  this  case  fl  is  decreasing 
on  (  —  1,0)  and  increa.sing  on  (0,1),  /z  is  increasing  on  (—1,0)  and  decreasing  on 
(0, 1),  fz  is  increasing  on  (—1, 1). 

C3:  Qif)  <  0  for  -1  <  e  <  0  and  0  <  ^  <  1:  In  this  case  fi  is  decreeising  on  (—1,1),  /2 
is  increasing  on  (  —  1, 1),  and  fz  is  increasing  on  (  —  1,0)  and  decreasing  on  (0, 1). 

C4:  Q{f)  >  0  for  — 1  <  ^  <  0  and  Q{f)  <  0  for  0  <  ^  <  1:  In  this  case  fi  is  increasing 
on  (  —  1,0)  and  decreasing  on  (0,1),  /z  is  decreasing  on  (  —  1,0)  and  increasing  on 
(0, 1),  fz  is  decreasing  on  (  —  1,1). 

C5:  Q(f)  =  0  for  — 1  <  ^  <  0  and/or  Qif)  =  0  for  0  <  ^  <  1.  In  this  case  /i,/2,/3  are 
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constant  on  the  region  where  Q{f)  =  0  and  /i,/2,/3  have  the  behavior  indicated 
in  Cases  1-4  where  Q{f)  >0  or  Q{f)  <  0. 

The  next  lemma  provides  L°°  and  total  variation  bounds  for  solutions  of  (Ve)- 

Lemma  1.2.  For  data  satisfying  (M),  the  functions  /i,  /2,  /s  are  positive  on  [-1,1]. 
Moreover,  /i,  /j,  /a  are  uniformly  bounded  from  above  and  below  by  positive  constants, 
and  of  uniformly  bounded  total  variation  on  [—1,1].  The  bounds  depend  on  the  data 
and  f~  but  not  on  e. 

Proof.  The  cases  C1-C5  are  ancilyzed  separately: 

Ci:  Here,  monotonicity  implies  0  <  /f  <  fi  <  <  f 2  <  ■  Since  Q  >  0  on 

[—1,1],  it  follows  that  /a  does  not  vanish  and 

0  <  <  iflhY^^  <  A  on  1-1,1], 

C2:  The  shapes  of  the  //s  in  this  case  dictate 

0  <  /,■  <  A  <  /,+  ,  0  <  min{/,-,/+  )  <  A  <  A(0) , 

A(0)  <  A  <  max{/f  ,/;*■)  on  |-1, 1] , 

Thus  /2  and  /a  are  positive.  Lemma  1.1  part  (ii)  states  that  Q  vanishes  at  ^  =  0, 
that  is  /i(0)/2(0)  =  /3(0),  hence  fi  is  also  positive.  Next,  use  the  balance  of  total 
mass  obtained  from  (1.8) 

-«A  +  A  +  4 A)l|iL,  +  (A  -  MlZt\  +  J\f,  +  A  +  4AK  =  o , 


15 


to  compute  the  total  mass 


J  (/i  + /z  +  =  2/j  +2/^+4(/3  +/^). 

Use  again  (1.8)  but  this  time  on  — 1  <  ^  <  0  to  arrive  at 

/2(0)  =  fih  +  /2  +  -  2/,-  -  if,-  +  /,(0)  <  C  , 

where  the  constant  C  depends  only  on  /*  but  is  independent  of  e.  Since  /i(0)  = 
/|(0)//2(0),  we  conclude  that  fi  and  /2  are  bounded  from  above  and  below  by 
positive  constants  depending  only  on  the  data. 

C3;  Here,  0  <  <  fi  <  fi,  0  <  <  f 2  ^  Lemma  1.1  part  (ii)  implies 

Qifi^))  =  0;  hence,  /3(0)  is  bounded  by  (/," and 

0  <  min{/3-,/+}  <f,<  (/r  . 

C.^:  The  same  argument  as  given  in  Case  2  applies  here. 

C5:  If  /  =  (/i,/2,/3)  is  constant  on  either  — l<^<0or0<^<l  the  monotonicity 
of  /  in  the  region  where  /  is  non-constant  immediately  yields  the  desired  bounds. 

The  above  arguments  show  the  stated  uniform  bounds.  Since  f\ ,  /2,  /a  may  have  at  most 
one  maximum  or  minimum  on  [—1, 1],  it  follows  that  /  =  (/i,/2,/3)  possesses  uniformly 
bounded  total  variation  as  well.  ■ 

Next,  we  examine  the  behavior  of  /  near  the  singular  points  with  the  goal  to  estimate 
the  modulii  of  continuity  of  the  components  fj.  To  this  end  we  use  representation  formulas 
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and  estimates  for  a  singular  linear  ordinary  differential  equation  that  are  derived  in  the 
Appendix.  Naturally,  it  is  not  expected  that  any  Holder  or  derivative  bounds  for  the  fj's 


will  be  independent  of  e.  To  differentiate  the  dependences,  for  the  rest  of  the  section  Cc 
will  stand  for  constants  that  depend  on  /*  and  e,  while  constants  that  depend  on  the  data 
but  are  £-independent  are  denoted  by  C. 

It  was  shown  in  Lemma  1.2  that 

0<5</>(O<C  {^1-1,1],  j  =  l,2,3;  (1.11) 

hence,  as  a  consequence  of  (0.16), 


M-i.i) 

i/;«)i  < 

^^[-l,l)\{0). 

(1.12) 


Relations  (1.12)  provide  derivative  bounds  away  from  the  singular  points  ^  =  —  1,  0  and  1 
respectively;  they  axe  improved  in  the  following  lemma. 


Lemma  1.3.  Let  /  =  (/i,/2,/3)  be  a  solution  of  (Ve)  with  e  >  0,  corresponding  to  data 
satisfying  (M).  There  is  pe,  1  <  p*  <  oo,  such  that  f[,  are  pj-integrable  on  [—1,1]. 

The  exponent  p*  depends  on  the  data  and  e,  it  is  increeising  as  e  decrezises,  and  there 
is  £o  =  £o(/^)  >  0  such  that  for  e  <  eo  the  exponent  pe  =  oo.  Moreover,  the  I^'-norms 
of  /{,  /j,  /3  are  bounded  by  constants  that  depend  on  and  e. 

Proof.  First,  we  work  with  the  component  /a  for  ^  >  0  in  a  neighborhood  of  0.  Recall 
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that  /a  eC[—l,  1]  satisfies 


-e/3  +  ^/3/3  =  ^/l/2  on  (0,1/2], 

2^  2£ 

/l(0)-/i(0)/2(0)=0 

With  the  purpose  to  use  the  results  in  the  Appendix,  note  that  p  —  fz  satisfies  an  equation 
of  the  form  (5.1)  on  (0, 1/2]  with  <t>  =  ^fz  and  h  =  Hypothesis  {H\)  requires 

W=  ^/3(0)>0, 


which  is  valid  in  our  case,  and 

—  /  ^  I  -^3(0  -  /3_(0)  I  (1  14) 

2^70  '  C  ' 

which  we  proceed  to  show. 

In  the  interval  (0, 1)  either  Q(f)  >  0,  or  Q{f)  <  0,  or  Q{f)  =  0.  For  concreteness,  the 
analysis  is  presented  for  the  case  Q  >  0.  (The  case  Q  <  0  is  treated  similarly  while  the  case 
Q  =  0  is  trivied).  Note  that  Q  >  0  implies  fz  increasing,  fz  >  (/i/z)^^^  and  fz  >  fz{0)  on 
(0, 1).  From  (1.13)  we  obtmn  the  identity 

Ml/2)  -  /,(0)  =  i  d( 

"  h ir  (  ~  +  VMOMO )  ii . 

which,  together  with  (1.11),  yields 


0  < 


MO- V MO f2(0 
< 


dC  <  eC 
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Using  the  fundamental  theorem  of  calculus,  (1.11)  and  (1.12),  we  deduce 

MO  -  MO) 


0< 


i: 


di 


_  MO-VMOMO  ,  r''WMOMO-VM0)M0) 


-I 


i: 


c 


dC 


^  MO- y/ flip f 2(0  0''^  (  nitOMtO  +  MtPMto  ^  \ 

Jo  C  Jo  \Jo  2  \/  fiitp  M^O  J 


VMiPMtP 

<  £  C  H —  . 

£ 

Thus  (1.14)  follows  and  Hypothesis  (/Ti)  is  satisfied. 

Observe  next  that  (5.6)  reads  in  the  case  of  (1.13): 

[MO  -  /3(0)]  +  /3(0)  ^^3(0)/2.  f  '  I^MO- MO  3-/3(0)/2.  gU,(£)-^(3) 

26  S 

=  [Ml/2)- MO)] 

^^M0)/2f  0'^  MOMO  -  /l(0)/2(0)  „-h(0)/2e  ^w(0-w(s)  ^,^5^ 

Jc  2s  s 


where 


u;(0  -w(s)  = 


MO  -  MO) 

2eC 


dC. 


The  two  first  terms  in  (1.15)  have  the  same  sign.  Thus  using  (1.15),  together  with  (1.11), 
(1.12)  and  (1.14),  we  find 

0  <  /3(O-/3(0)  <  j  '  5-A((I)/2.J3 

^  f  if  MO)  ^2s  J 

-\c.{(H-lKi)  if/3(0)  =  2£ 

Equation  (1.13)  implies  with  the  help  of  (1.11),  (1.12)  and  (1.16)  that  for  ^  e  [0, 1/2] 


0  <  1/^(0!  = 


2^[(/l  -  fim  -  (hh  -  /.(oj/jco))] 


(l  +  kl^-')  if/3(0),i2£ 


(1  + 


if/3(0)  =  2£ 


(1.17) 
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This  improves  (1.12)3  for  ^  >  0  in  a  neighborhood  of  0.  Since  f3{-()  satisfies  an  equation 
of  the  form  (1.13)  as  well,  the  same  analysis  shows  that  (1.17)  in  fact  holds  in  a  full 
neighborhood  of  0  and  ultimately,  upon  combining  with  (1.12)3,  for  ^e[— 1,1].  From 
(1.17)  we  estimate  the  Z^-norm  of  /j  on  the  interval  [—1, 1] 

(i)  If/3(0)>2£  then||/'lU«  <C.. 

(ii)  If  /3(0)  <  2£  then  for  pe  <  1/(1  -  we  have  ||/3||z,p.  <  Cj. 

Next  we  turn  to  the  component  fi .  If  we  set  p(0  =  /i(l  -0  then  p  satisfies  an  equation 
of  the  form  (5.1) 

-^p'(e)  +  i/2(l-^)p(^)  =  i/|(l-^)  on  (0,1/21, 

P<0)  =  /,+  =  fUl)/Ml) 

Hypotheses  {H\  —  H2)  are  satisfied  with  /?  =  7  =  1  and  (5.12),  (5.16)  together  with  (1.11), 
(1.12)  yield 

l/!(el  <  I  tor  Ce|-l,l]  (1.18) 

lC.(l  +  Ml-OI)  if/?=£ 

In  turn,  (1.18)  implies  that  ||/|||i,j>c  <  Ce  for  =  00  in  case  >  £,  and  for  any 
Pe  <  1/(1  -  ^)  in  case  <  e. 

Finally,  a  similar  argument,  working  with  the  function  p(^)  =  /2(— 1  +  0?  shows  that 
1/2(01  <(  for  ^e[-l,l]  (1.19) 

I  Ce(l  +  |/n(l  +  OI)  if/i"=£ 

and  that  ||/2||i,p.  <  Ce  for  p^  =  00  in  case  /j~  >  £,  and  for  any  pe  <  1/(1  —  ^)  in  case 

fr 
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We  emphasize  that  for  e  <  min{/f ,  /3(0)/2}  the  derivatives  /j(^)  axe  uniformly 

bounded  on  [—1,1].  Since  /3(0)  is  bounded  from  below  by  (1.11),  there  is  a  threshold 
eo  >  0  (depending  only  on  /^)  so  that  for  e  <  Sq  the  norms  ||/j||i,~  <  Cg,  j  =  1,  2,  3. 
Naturally,  the  bounds  blow  up  as  £  — »  0.  i 

A  useful  consequence  of  the  derivative  bounds,  the  Sobolev-type  inequality 

l/>(«-/i(C)l  =  I  s  ll/ill'.’- 

and  Lemma  1.2  is  the  following  corollaury: 

Corollary  1.4.  The  functions  /i,  /2,  /3  axe  Holder  continuous  with  exponent  Og  =  1  -  ^ 
(ag  =  1  if  pg  =  oo).  Moreover,  the  Og-Holder  norms  of  /i,  /2,  fz  are  bounded  by  constants 
depending  on  and  £. 

We  close  this  section  with  a  compactness  lemma. 

Lemma  1.5.  Let  £  >  0  be  fixed  and  let  {/*'}  =  { (/f ,  /2  ,/3*)}  be  a  sequence  of  solutions 
to  the  nonlinear  boundary  value  problem  (“Pg)  taking  on  Maxwellian  data  at  ^  =  ±1. 
Assume  the  data  are  uniformly  bounded  from  above  and  below  independently  of  k,  i.e. 
there  exist  constants  0  <  6,  M  <  oo  so  that  6  <  f  ^  ^  M .  Then  there  exist  a 

subsequence  {f^  }  of  {/*}  and  a  continuous  function  /,  such  that  f*‘  f  uniformly  on 
[-1,1],  and  /  satisfies  (Pg)  and  admits  positive  Maxwellian  boundary  data. 

Proof.  From  Corollary  1.4  we  see  that  for  each  fixed  j  the  sequences  {fj}k>i  li®  i^^ 
bounded  sets  of  C“[— 1, 1].  Since  on  [—1, 1]  the  embedding  C®  — »  (^  <  o)  is  compact. 
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there  is  a  subsequence  {f'‘  }  and  functiono  fj  eC"  such  that  fj  — ♦  fj  in  .  In  particular, 

and  since  the  data  of  the  sequence  correspond  to  positive  Maxwellians,  the 
boundary  values  /(±1)  =  will  be  positive  Maxwellians  as  well.  Passing  to  the  limit 
k'  oo  \n  (1.1-1. 3)  we  see  that  /  is  a  solution  of  {Ve)  (in  the  sense  of  the  definition)  zind 
thus  enjoys  all  the  additional  regularity  properties  that  were  proved  in  this  section.  | 

§2.  The  fluid  dynamic  limit. 

Let  be  fixed  positive  Maxwellians  and  let  {f‘}e>o  be  a  family  of  solutions  of  (Ve) 
admitting  /*  as  boundary  data.  The  members  of  the  family  are  extended  to  the  whole  real 
line,  by  setting  /®(^)  =  f~  for  ^  — 1  and  f‘{^)  =  /"*■  for  ^  >  1.  The  extended  functions 
are  again  denoted  by  f‘ .  In  this  section  we  show  that  the  family  of  the  extended  functions 
possos.ses  a  subsequence  {/*’’}  with  £„  — ♦  0  which  converges  pointwise  to  a  solution  /  of 
the  fluid  dynamic  Riemann  problem. 

Theorem  2.1.  Let  {/'}e>o  be  a  family  of  solutions  of  (Ve)  corresponding  to  data 
satisfying  (M).  There  exists  a  subsequence  {/'"}  with  £„  — »  0  and  a  positive,  bounded 
function  f  of  bounded  variation  such  that  /'•*  — »  /  pointwise  on  the  reals.  The  function  f 
is  a  local  Maxwellian,  that  is  ff  =  /1/2  for  a.e.  and  satisfies  the  balance  of  mass  and 
momentum  equations 

-f^(/l+/2+4(/,/2)''n  +  ^(/.  -h)  =  0.  (2.1) 

+  +/2)  =  0,  (2.2) 

in  the  sense  of  distributions  and  in  the  sense  of  measures. 
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Proof.  Since  Q{f^)  =  0  the  functions  /*  satisfy  the  equations  (1.1-1. 3)  for  any  ^i,^2  in 
the  reals.  Lemma  1.2  in  conjunction  with  Helly’s  theorem  implies  that  {f‘}  possesses  a 
subsequence  {/*"}  with  £„  — »  0  which  converges  pointwise  on  [—1,1]  to  a  function  /  = 
(/i,/2,/3)  with  positive,  bounded  components  of  bounded  variation.  The  functions 
admit  constant  values  /*  outside  [—1,  Ij;  therefore,  first  /*"  —*  f  pointwise  on  (—00,00), 
and  second  /(^)  =  f~  for  ^  —  1  and  f(^)  =  /"*’  for  ^  >  1- 

Let  now  ^1,^2  e  [— 1,  !]•  Equation  (1.10)  for  the  functions  /*",  with  the  help  of  the 
inequality 

(a  -  b)iCna  -  inb)  >  4(a‘/^  -  ,  a,  6  >  0, 


and  the  uniform  bounds  (1.11),  gives 

r?2 


C  -  /f" )(fn  ft-^-  fe/,'" /'" ) 

hi 

<  -4  /I-)''")"  <'«<». 

hi 

Passing  to  the  limit  £„  — >  0,  we  find  from  the  dominated  convergence  theorem 

hi 


(2.3) 


and  thus  /j  =  (/i/z)’^^  a.e.  in  [—1,1]. 

Let  ip{0  be  a  function  with  compact  support  in  (—00,00).  Equation  (1.5)  gives 


/"(/r-  +/2'"  +4/'")(^9(o)'-(/r  -f2’')^'iOd^  =  o- 

J —00 

Letting  £„  — »  0  and  using  once  again  the  dominated  convergence  theorem,  we  obtain 

/°°(/i  +  A  +4(/,/2)'/")  («,»{{))' -(/■  =  0.  (2.4) 

J  —00 
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Thus  (2.1)  holds  in  the  sense  of  distributions  and,  because  fi,  /z  are  of  bounded  variation, 
it  also  holds  in  the  sense  of  measures.  Parsing  to  the  limit  t„  — ♦  0  in  (1.6),  we  show  (2.2) 
and  complete  the  proof.  | 

With  fi,  f 2  as  above,  define 


Fi{x,t)  =  fx{j),  F2(i,t)  = /2(j),  (x,t)e(-oo,oo)  X  (0,oo)  (2.5) 

Clearly  limt_o  Fj{x,t)  =  f~  for  x  <  0,  fj'  for  x  >  0,  j  =  1,  2.  Further,  a  solution  (Fi,  F2) 
of  the  form  (2.5)  is  a  weak  solution  of  (0.6)  on  (—00,00)  x  (0,oo)  if  and  only  if  (/i,  /2)  is 
a  weak  solution  of  (2. 1-2.2)  on  (-  xd,oo). 

The  equivalence  follows  from  a  transformation  among  test  functions  (Dafermos  [7]). 
Indeed,  let  rl^{x,t)  be  a  C°°  function  with  compact  support  on  (—00,00)  x  (0,oo)  and 


define' 


^00 

S?(0  =  / 

Jo 


(2.6) 


The  resulting  function  ^  e  and  has  compact  support  on  (— cx5,  00).  Conversely,  any  test 
function  (p  may  be  represented  in  the  form  (2.6)  by  choosing  V’  =  with  a{t) 

a  fixed  C°°  function  compactly  supported  in  (0, 00)  and  satisfying  a{t)tdt  =  1.  For 
solutions  of  the  type  (2.5),  the  weak  formulation  for  the  first  (sa> )  equation  in  (0.6)  takes 
the  form 

fOO  /*oo 

/  /  (Fi  -HF2-t-4(F,F2)‘/2)0,(x,O  +  (Fi  -F2)0.(x,f)dxdf 

Jo  J —00 

/oo  roo  roo 

(/l  +/2+4(/,/2)‘/^)(  /  M^t,t)tdt)+(fi-f2)(  /  M^,t)tdt)d^ 

=  /”(/.  +  /2  +  4(/,/j)‘e^)  (  -  «9(0)'  +  (/i  -  h)v\0di 

J  — OO 
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eind  the  equivalence  follows  from  the  chain  of  identities. 


§3.  Existence  theory  for  the  nonlinear  boundary  value  problem  {Ve) 

The  scope  of  this  section  is  to  develop  an  existence  theory  for  the  nonlinear  boundary 
value  problem  (Te),  consisting  of  (0.16)  and  (0.18),  for  £  >  0  fixed.  Throughout,  the 
boundary  data  /*  are  assumed  to  be  positive  Maxwellian  states  satisfying  (Af). 

First,  some  convenient  notation  is  introduced.  Let 

K/)=(-/2,-/i,2/3)  (3.1) 

and  define  the  matrices 

— /2  — /i  2/3 

B(/)=  -h  -/i  2/3  (3.2) 

.\h  kh  -h. 

0  O' 

MO=  0  -(?  +  l)  0  ,  -1<{<1.  (3.3) 

0  0  -i 

Then  Q{f)  =  ^b{f)  •  /,  and  (0.16)  may  be  written  in  the  form 

■  (3.4) 

Next,  let  F  =  {Fi,  F2,  F3)  be  a  given  solution  of  (Ve)  defined  on  [—1,1]  and  taking 
boundary  values  F(±l)  =  F^  positive  Maxwellian  states.  The  analysis  in  Section  1  implies 
that  the  functions  Fj  are  Holder  continuous  with  exponent  Oei  moreover,  for  £  <  £o{F^) 
they  are  Lipschitz  continuous.  Any  nearby  solution  /  can  be  written  as 

f  =  F  +  p, 

25 


with  p  —  (pi,P2,P3)  a  perturbation.  Since 


Qif)  =  Q(F  +  p)  =  Q(F)  +i(F)-p  +  Q(p)  ,  (3.5) 

•t. 

p  is  a  continuous  function  on  [—1, 1]  that  satisfies  the  equations 

A(Op’  =  \b(F(,())p  +  ^B(p)p  ,  (3.6) 

for  ^  €  (  —  1,0)  and  ^  €  (0, 1),  and  the  boundary  conditions 

p(±l)=p±:=/±-F±.  (3.7) 

Since  /*  and  F*  are  Maxwellian  states,  (3.5)  imposes  the  restriction 

6(F±)-p±  +  Q(p±)  =  0  (3.S) 

on  the  data  and  that  in  turn  implies 

B(^'*)p*  +  |b(p*)p±=0.  (3.9) 


Our  strategy  for  proving  existence  is  to  use  a  continuation  argument  on  the  set  of  data 
for  which  (Ve)  admits  a  solution.  Any  constant  Maxwellian  is  a  trivial  solution  of 
(Vs),  corresponding  to  data  with  and  can  serve  as  a  point  of  departure  for 

the  continuation  argument.  The  key  ingredients  are:  Lemma  1.5  and  showing  that  for 
prescribed  F(^)  and  with  —  F*  sufficiently  small  a  solution  of  (3. 6-3. 7)  exists.  This 
second  objective  is  pursued  here. 
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We  seek  for  solutions  of  (3. 6-3. 7)  as  fixed  points  of  the  map  that  carries  a  continuous 
function  P  =  (Pi,P2,  P3)  with  boundary  values  P(±l)  =  p^,  to  the  solution p  =  {p\,P2iPz) 
of  the  boundary  value  problem 


mp'  =  + ^B(P((mo 


(3.10) 


p(±l)  =  p^  . 

The  study  of  this  map  is  based  on  properties  of  linear  boundary  value  problems  near  the 
regular  singular  points  ^  =  0,  ±1,  that  are  developed  in  the  following  subsection. 


§3.1.  The  linearized  boundary  value  problem 
Consider  the  singular,  linear,  boundary  value  problem 

X(0p'  =  iB(r(0)p  +  iG({)  ,  -1<{<1, 

p(±l)  =  p*  . 


(3.11) 

(3.12) 


The  matrix  B(F(^))  is  defined  in  (3.2)  and  has  Holder  continuous  entries  on  [—1,1]  with 
exponent  Og,  while  G  =  (<71,  <72 1  <73)  is  a  continuous  vector  valued  function  on  [—1,1].  The 
boundary  data  p^  are  assumed  to  satisfy 

6(P'’’)  •p'*'  -byi(-f-l)  =  0 


(3.13) 


b(F  )-p  •hff2(-l)  =  0  . 

Note  that  the  boundary  value  problem  (3.10)  fits  into  this  framework  with  G  =  ^B(P)P, 
and  that  relations  (3.13)  refiect  the  restrictions  imposed  by  (3.8). 

Our  goal  is  to  study  the  solvability  of  (3.11-3.12)  subject  to  (3.13),  and  to  establish 
estiii'  "tes  for  the  solution  p  in  terms  of  the  data  G  zind  p^.  For  the  estimates  we  use  the 
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sup-norm  ||  •  ||  for  scalar  or  vector-valued  functions  (depending  on  the  context).  The  results 
axe  summzirized  in  the  following  Fredholm-alternative  type  of  theorem. 


sup  |p(OI<C'(  sup  |(?(0I  +  +  b  !)•  (3-15) 

Proof.  The  proof  is  long  and  will  be  established  in  four  steps.  The  main  obstacle  is  the 
presence  of  singularities  at  ^  =  —  1,  0  and  1. 

The  first  two  steps  are  preparatory  in  nature  and  concern  the  behavior  of  solutions  in 
the  neighborhood  of  one  singular  point.  It  is  expedient  to  consider  the  auxilliary  system 

-«p;  =  -i'Pcejpj  +  ;<«)  •  (p.,P2) + ,  (3.17) 

taken  on  the  interval  [0,  |]  under  the  list  of  hypotheses:  o  is  a  nonsingular  continuous 
matrix  on  [0,  |],  0  a  continuous  matrix,  7  =  (71,72),  ^  =  (^1,^2)  and  H  =  (/ii,/i2) 
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continuous  vector-valued  functions,  /i  is  a  continuous  scalar-valued  function,  and  9?  >  0  is 
a  Holder  continuous,  positive,  scalar-valued  function. 


Step  1.  Behavior  of  solutions  near  a  singular  point. 

It  is  well  known  that  for  a  linear  homogeneous  system  of  nonsingular  ordinary  differential 
equations  the  solution  operator  is  an  isomorphism  in  the  Euclidean  space.  We  proceed  to 
study  the  analog  of  this  result  for  a  singular  system  with  the  structure  of  (3.16-3.17). 

In  preparation,  consider  (3.16-3.17)  supplemented  with  the  initial-boundary  conditions 

Pi(0)  =  Pio  ,  P2(0)  =  P20  ,  P3(a)  =  P3a  ,  (3.18) 

where  the  last  condition  is  applied  at  some  intermediate  point  a  €  [0,  ^].  We  show: 

Theorem  3.2.  It  a  is  sufficiently  small,  there  exists  a  unique  continuous  solution  of 
(3.16-3.17)  that  admits  the  data  (3.18).  It  satisfies 

<P(0)P3(0)  -  <5(0)  •  (pi(0),p2(0))  =  h{0) .  (3.19) 

Proof.  We  apply  the  contraction  mapping  theorem  to  an  equivalent  integral  formulation 
of  (3.16-3.18). 

Let  ^(0  t>e  a  fundamental  matrix  solution  for 

(::)'= 

with  ’J'(O)  =  I.  The  variation  of  constants  formula  applied  to  (3.16)  implies 

\j‘2j  £  7o  (320) 

+  -<6(0  f  <i-\s)a-'(s)H(3)ds  . 

^  Jo 
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Thr'  fot  iuulft  (5,fl)  in  thn  Appninlix,  ii|»plir^«l  in  (!M7),  yitiMa 


wln^rn  in  ihia  I'ttaf 


Ji 


=  iil  (. 


I'hn  inlrmnl  r**piniiona  (ll.lJt)  a.'Jl)  pmvidr’  nti  n»|nivnlnni  fninnilniinn  for  ihn  pjohlr^in  >’t>n 
aiaiinH,  of  {3. 13  3,  IS), 

N*’Xi,  aoi 

ViO-vPisM  /  vl'  '(a)o-'(a)/f(a)./a 

Ji\ 

u 

tnnl  tlt'ilnf’  iliti  innpa  T  nn«l  sV,  on  apno^a  t>f  *'oniinno»»a  fnin'iiona,  na  followa; 

(i)  Tin*  mnj)  T  cnrrit^a  ilit^  pair  of  ooniinttona  fmu’iiona  {A,  A)  Pi  Ktv*^n  hy 

(I 

^  Ji  «= 

If  (/’i,/’ji),  {P\,Pi)  tw«>  paira  of  ooniinnona  fnin'iiona  and  p:\  —  T(P\^f\),  - 

T{P\,  Pi)  ihair  r»’apt't’tiva  iinagna,  ihan 

U'<  <(»)■  ((/',  ^P, )(..), (/'j  -riX,.)),, 

anti  a  tlat ivaiion  in  ilm  apirii  t\f  Lamina  ft,  I  imliaaiaa  ihai  tha  map  T  aaiiaHaa  ilia  aaiimaia 


p.  -  P:,ll  ^  a(i|/-,  -  P,||  +  II  n  -  P,|| , 


3(1 


where  ||  •  ||  stands  for  the  C[0,  a]  sup-norm  and  the  constant  C  is  independent  of  a. 

(ii)  The  map  S  carries  the  continuous  function  P3  to  the  pair  (pi,P2)  defined  by 

(!')  =  f  Pz(s)<i-\s)a-'(s)~,(s)ds  +  . 

\P2/  \P20/  £  Jo  £ 

If  now  P3  and  P3  are  two  continuous  functions  and  (pi,P2)  =  S{P3),  (Pi,P2)  =  ‘^(Ps) 
their  respective  images,  then 

fP>  ~Pj\  =  f  \p,-7,)(s)i-'(s)a-'{s)y{s)d3  , 

\P2  -  P2J  £  Jo 

and  it  follows  easily  that  the  map  S  satisfies  the  estimate 

IIPi  -  Pill  +  IIP2  -P2I!  <  Ca||P3  -  P3II  (3.24) 

with  the  constant  C  independent  of  a. 

Combining  (3.23)  with  (3.24)  we  see  that  the  composite  map 

S  oT  :  C[0,  a]  x  C[0,  a]  — +  C[0,  a]  x  C[0,  a] 

is  a  contraction,  provided  that  a  is  sufficiently  small.  Therefore,  S  oT  admits  a  unique 
fixed  point  which  is  a  solution  of  (3.20-3.21).  Finally,  (3.19)  follows  from  relation  (5.7)  in 
the  Appendix.  | 

Consider  now  the  homogeneous  system  of  (3.16-3.17)  with  P  =  0,  h  =  0.  Let  4>(2) 
and  <^(3)  be  the  solutions  corresponding  to  the  choices  of  data  (pio,P20  5P3o)  equal  to 
(1,0,0),  (0, 1,0)  and  (0,0, 1)  respectively.  Then  ^(2),  ^(3)  form  a  fundamental  set  of 
solutions  for  the  homogeneous  system  (3.16-3.17).  They  also  enjoy  the  properties: 


31 


(a)  For  ^  ^  0  the  vectors  <l>(2){0^  <^(3)(0  ^xe  linearly  independent. 

(b)  For  ^  =  0  the  vectors  <^(i)(0),  <^(2)(0),  ^(3)(0)  are  linearly  dependent  and  span  the 
two-dimensional  space  defined  by  (3.19)  for  h  ~  0. 

The  general  form  of  a  fundamental  matrix  for  the  homogeneous  (3.16-3.17)  is  given  by 

[</»(!)  , 

where  each  <l)j  is  a  column  vector  and  5  is  a  nonsingular  constant  matrix.  From  the  form 
of  the  general  fundamental  matrix,  it  follows  that  any  fundamental  set  enjoys  properties 
(a)  and  (b). 

Step  2.  An  a-priori  estimate 

Let  now  p  =  (pi  ,P2,P3)  be  a  solution  of  (3.16-3.17)  on  (0, 1]  (not  necessarily  defined  at 
^  =  0),  and  introduce  the  notation  ry  =  j  =  1,2,3.  We  show 

Lemma  3.3.  The  solution  p  can  be  extended  to  a  continuous  function  on  [0,  |]  that 
satisfies  (3.19)  and  the  estimate 

H<c(^|r,|  +  ||^r||  +  ||Ml),  (3.25) 

J=1 

where  ||  •  ||  stands  for  the  (vector  or  scalar)  sup-norm  on  [0,  |]. 

Proof.  The  variation  of  constants  formula  yields  from  (3.16) 

(^‘)  =  «(e'f'(i)(|;^)  -  -  ic/(e ,  (3.26) 

where 

/•1/2 

U{0  =  ^(0  J  '^-\s)a-^(s)H{s)ds  .  (3.27) 
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From  (3.26)  we  obtain 


|P.(0I  +  |P2{«)I  <  c(|r,H-  Ir^l  +  y  '  |p3(^)|*  +  lU(Ol) 


Next,  (3.21)  gives  for  a  =  | 


P3  = 


where 


Using  the  estimate 


,«i;(^)-u>(a)  _  7  // 


<  C*  Jo 


(3.28) 


(3.29) 


(3.30) 


the  inequality 

^iv>(o)  (3  31) 

V’(O) 

jmd  (3.28),  we  deduce 

|P3{«I  <  C(|r,|  +  \u(()\)+C(i'^°'  j\\pi(s)i  +  |p2{s)|)s-J'l'’)-‘ds 

<c(E|r,|  +  K«|)  + 

j=l 

+  C{i»<”>  J’  (/^P3(C)l<^C)»-•''‘"’■‘<i» 

pi 

^  yi 

<C  (E  hi  +  lie'll  +  ll'‘ll) +c  /  MOMC- 

>=1 
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Gronwall’s  inequality  then  implies 

3 

|P3(f)l  <  C(^  |ry|  +  llt^ll  +  Hull)  ,  0  <  f  <  i  .  (3.32) 

Therefore  is  bounded  and  thus  pi,  p2  are  Lipschitz.  Lemma  5.1  then  implies  that  P3(0 
has  a  limit  as  ^  0+  and  the  limiting  functic.  satisfies  (3.19).  As  a  matter  of  fact, 

Lemma  5.2  implies  that  pa  is  Holder  continuous  on  [0,  |]  On  account  of  (3.28)  and  (3.32), 
we  obtain 

3 

\M0\  +  |P2(0I  <  c  (E  hl  +  ||t^ll  +  l|u||). 

>=1 

Finally,  (3.27),  (3.30)  and  (3.31)  imply 

ii?7ii  <  cm ,  ii«ii  <  cii/i|i 

and  the  proof  is  complete.  § 

Step  3.  The  singular  system  (3.11) 

From  now  on  we  focus  on  the  system  (3.11).  First  we  explore  the  relation  between  (3.11) 
and  the  auxilliary  system  (3.16-3.17)  studied  in  Steps  1  emd  2.  Clearly,  the  latter  captures 
the  local  structure  of  the  former  in  a  positive  neighborhood  of  0  and  thus  describes  the 
behavior  of  solutions  to  (3.11)  £is  ^  ♦  0-I-.  More  is  true:  (3.16-3.17)  captures  the  local 

structure  of  (3.11)  in  the  neighborhood  of  each  of  the  singular  points,  not  only  as  ^  »  0-f 
but  also  as  ^  0—,  as  ^  »  1—  and  as  ^  — 1-t-.  To  see  that,  consider  a  solution 

(Pi(0»  P2(0>  P3(0)  (3-11)  and  suppose  we  are  interested  in  the  behavior  as  ^  1—. 

Performing  the  change  of  variables  ^  =  1  —  C,  we  deduce  that  the  triplet 

(?l(C).  «(C).  ?3«))  =  (P2(l  -  C).  P3(l  -  C).  Pl(l  -  0) 
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satisfies  for  0  <  C  ^  2  ^  system  of  the  form  (3.16)-3.17)  with  <^(C)  =  ^2(1  —  0  >  0  ^(C) 

a  nonsingular  matrix  on  [0,  |].  Lemma  3.3  implies  that  p  can  be  extended  to  a  continuous 
function  as  ^  1— ,  that  the  limiting  value  satisfies 

6(F(+l)).p(+l)  +  ^i(+l)  =  0, 

and  that  p  can  be  estimated  by 

3  .. 

sup  Ip(0l  <  bi(+«)l  +  sup  |G(0l)  • 

Similar  statements  follow  from  analyzing  the  behavior  as  ^  +  0—  and  as  ^  »  —1+-  They 
are  summarized  below: 

Corollary  3.4.  Let  p  be  a  solution  of  (3.11)  on  (—1,0)  U  (0, 1).  Then  p  can  be  extended 
to  a  function  that  is  Holder  continuous  on  each  of  the  interv'als  [—1,0]  and  [0, 1]  (though 
not  necessarily  continuous  at  ^  =  0)  and  which  satisfies: 

(i)  At  the  singular  points  the  conditions 

6(F(+l))-p(+l)  +  pi(+l)  =  0, 

6(F(-1))-p(-1)  +  J2(-1)  =  0,  (3.33) 

i<F(0))-p(0±)+j,(0)  =  0. 

(ii)  The  estimates 

I|p||±<c(^|r,±|  +  ||G|U),  (3.34) 

7=1 

where  Tj±  =  Pj{±\)  ,  j  =  1,2,3  ,  ||  •  ||+  stands  for  the  sup-norm  on  [0, 1]  and  ||  •  ||-  for 
the  sup-norm  on  [—1,0]. 


Consider  now  the  homogeneous  system  (3.11)  with  G  =  0.  Choose  a  suffiently  small  and 
assign  data  for  (pi(0),  P2(0),  P3 (a))  equal  to  (1,0,0),  (0, 1, 0)  and  (0, 0, 1).  From  Theorem 
3.2,  we  obtain  a  fundamental  set  of  solutions  on  [0, 1].  On  account  of  (3.33) 

the  functions  take  boundary  values  satisfying 

i<F+)  •  .^+,(+1)  =  0  .  6(f(0))  •  <^,(0+)  =  0  .  (3.35) 

Similarly,  by  assigning  data  for  (pi(0),  P2(0),  P3(— a))  equal  to  (1,0,0),  (0,1,0)  and 
(0,0,1)  we  construct  a  fundamental  set  of  solutions  ^(3)’  respectively,  on 

[—1,0].  The  boundary  values  of  these  functions  now  satisfy 


b(F.)  •  <^-^(-1)  =  0  ,  6(F(0))  •  fli(-)(0-)  =  0  .  (3.36) 

By  their  construction  and  (3.35),  (3.36), 

<,)(0+)  =  '!‘r,)(0-) ,  <2,(0+)  =  <2,(0-) ,  <3|(0+)  =  «1(-3|(0-)  =  0  .  (3.37) 

Next,  we  turn  to  the  inhomogeneous  system  (3.11)  and  construct  a  particular  solution 
on  [0, 1]  by  solving  (3.11)  subject  to  the  initied  condition  0'*’(+^)  =  0.  Note  that  the 
resulting  solution  is  bounded,  continuous  and  takes  values  at  ^  =  0+,  +1  satisfying 
(3.33).  Similarly,  we  construct  a  particular  solution  on  [—1,0]  which  satisfies  = 

0  and  the  conditions  (3.33)  at  ^  =  —1,0—.  Corollary  3.4  and  the  construction  also  dictate 


||#+|4  <  C||G||+  ,  IIV'II-  <  C||G||-  . 


(3.3S) 


The  general  solution  of  the  inhomogeneous  system  (3.11)  is  given  by 

f  +  02^6- )(^)  +  ,  -1  <  ^  <  0 


P{0=< 


+  63<i+3,(0  +  rk^iO  ,  0  <  ^  <  1 


(3.39) 
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with  ai,  02,  03,  bi,  63,  63  2irbitrary  constzints. 

Step  4-  singular  boundary -value  problem  (3.11-S.12) 

The  boundary  conditions  (3.12)  imply 

+  «2<^("2)(-1)  +  ^34>i3)(-'^)  =P--  , 

^l<^j'’l)(  +  l)  +  ^2<?i^)(  +  l)  +  ^3<^^)(  +  l)  =  P+  -  l/’'^(  +  l)  . 

On  the  other  hand,  the  requirement  of  continuity  p(0+)  =  p(0— )  =  po  dictates 

+  ^2<f>^2)(^-)  +  «3</^('^)(0-)  =  Po  -  t^"(0-)  , 

+  )  +  ^2</»^)(0+)  +  &3<^(^)(0+)  =  Po  —  ^^'*’(0+) . 

In  view  of  (3.37),  it  is  equivalent  to  solve  the  inhomogeneous  algebraic  system 
ai^(“i)(-l) +a2<^(l)(-l) +  a3<^(l)(-l)  =P-  -rp~i-l)  , 

(oi  -  6i)<?(~i)(0-)  +  (02  -  62)<;^(1)(0-)  =  i/>''‘(0+)  -  V’~(0-)  ,  (3.40) 

+  ^2<^^)(  +  l)  +  ^3<^^)(  +  l)  =  p+  —  0'^(  +  l)  . 

This  is  a  system  of  nine  equations  in  six  unknowns.  However,  because  of  (3.33)  applied 

to  the  functions  tjj~,  the  relations  (3.35-3.36)  and  the  compatibility  conditions  (3.13), 

at  most  six  equations  are  linearly  independent.  Solvability  of  (3.40)  for  any  choice  of  the 

data  p±  is  equivalent  to  the  homogeneous  algebrciic  system  having  only  the  trivial  solution, 

which  is  in  turn  equiveilent  to  the  homogeneous  problem  (3.14)  possessing  only  the  trivial 

solution  p  =  0  in  the  class  of  continuous  on  [—1, 1]  functions.  This  completes  the  proof  of 

part  (i)  of  the  theorem. 

To  prove  part  (ii)  observe  that,  if  the  homogeneous  (3.40)  has  only  the  trivial  solution, 

then  the  solution  (oi,02,03,bi,b2,b3)  of  the  inhomogeneous  (3.40)  satisfies  the  estimate 
3 

Y,  2  Clip-  -  +  IV’'"(0+)  -  ^-(0-)|  +  |p+  -  .  (3.41) 

>=1 
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Moreover,  evaluating  (3.39)  at  ^  =  ±^,  we  deduce 

3  3 

^  kj+l  +  k>-|  <  kjl  +  |6j|  .  (3.42) 

>=1  i=i 

Finally,  combining  (3.34)  with  (3.41),  (3.42)  and  (3. 38),  we  arrive  at  (3.15).  | 

§3.2.  The  nonlinear  boundary-value  problem 

Next,  we  study  the  solvability  of  (3. 6-3. 7).  Recall  that  t  >  0  is  fixed,  F{^)  is  a  given 
solution  of  (Ve),  while  the  data  p*  satisfy  the  restriction  (3.S).  We  show: 

Theorem  3.5.  If  the  only  continuous  solution  of  (2.1  A)  is  the  trivial  solution  p  =  0,  there 
exists  a  positive  constant  r  such  that  for  any  data  satisfying  (3.8)  and  |p'''|  -f-  |p~|  <  r  the 
boundary -value  problem  (3. 6-3. 7)  has  a  solution. 

Proof.  Since  our  proof  will  rely  on  the  contraction  mapping  theorem  first  define  the  closed, 
bounded  subset  of  the  continuous  functions 

A={Pe  C[-l,  1]  :  P(±l)  =  ,  ||P||  <  m} 

where  |j  •  ||  stands  for  the  sup-norm  and  m  is  a  positive  constant.  Consider  the  map  T 
that  carries  P  €  to  the  solution  p  of  the  boundary- value  problem  (3.10).  In  view  of  the 
Theorem  3.1  and  the  compatibility  conditions  (3.8)  the  map  T  is  well  defined. 

Our  goal  is  to  show  that  upon  choosing  m  and  r  sufficiently  small  the  map  T  •.  A —*  A 
and  is  a  contraction.  The  resulting  fixed  point  p  ^  A  will  then  be  the  claimed  solution. 

To  accomplish  that,  observe  that  (3.15)  implies 

l|p||<C(|p-^|  +  |p-|  +  ||iB(i’)F||) 

<  A'ldp-^l  +  Ip'I  +  lien’)  <  A-,(r  +  m’)  . 
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On  the  other  hand  if  P,  P  e  A  and  p  =  T{P),  p  =  P(P)  then  using  again  (3.15)  we 
conclude 

l|p-p||<C||iB(P)P-iB(P)P|| 

<  A-j(||P||  +  ||P||)||P  -?||  <  2K,m\\P  -  P||  . 

Choosing  m  <  min{y^,  2;^}  then  r  <  we  see  that  ||p||  <  m  and  that  ^  :  A  A 

is  a  contraction.  ■ 

§3.3.  Nonexistence  of  eigensolutions 

In  order  to  apply  Theorem  3.5,  we  need  to  check  that  the  homogeneous  boundary- value 
problem  (3.14): 

-K-i)p'i  = 

-«+l)P2  =  b(P(«))-P  (3.43) 

-(pz  =  -~mo)-P 

with  boundary  conditions  p(±l)  =  0,  does  not  admit  any  nontrivial  continuous  eigen- 
solutions.  We  have  been  unable  to  do  that  for  a  general  solution  P(^)  of  (Pc).  (Note 
however  that  such  a  result  is  true  in  the  case  of  the  Carleman  system;  see  the  next  Sec¬ 
tion.)  Nevertheless,  it  is  possible  to  rule  out  eigensolutions  for  certain  classes  of 
namely: 

Lemma  3.6.  In  either  of  the  cases: 

(a)  P(^)  is  a  constant  Maxwellian  state,  or 

(b)  g(P(0)>0  fore  €[-1,1], 

the  homogeneous  problem  (3.14)  possesses  only  the  trivial  solution  p  =  0  in  the  class  of 
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continuous  functions. 


Proof.  The  proof  is  based  on  an  energy  identity  for  solutions  of  (3.43).  Multiplying  the 
first  equation  in  (3.43)  by  F2P1,  the  second  by  F1P2,  the  third  by  4F3P3  and  adding  the 
resulting  equations,  we  obtain 


-(f  -  m  (5?;)'  -  (?  +  l)f.  (|pi)'  -  4«/'3(|p§)'  =  -i(6(i")  •  P)“ 


Integrating  by  parts  and  using  (0.16),  we  arrive  at 

+  fijpl  +  2F3Pi)'  +  (f^P?  -  •F'i5Pi)'+ 


+ 


(3.44) 


We  integrate  (3.44)  over  [—1,1]  and  use  the  continuity  of  p  at  ^  =  0  and  the  boundary 
conditions  p(±l)  =  0  to  conclude 

J  (^o^2Pi  +  ;^Fip2  4- 

3  /' (fTT)"?  +  3  =  « • 

Since  F  is  a  solution  of  Vc,  the  functions  Fj  are  strictly  positive.  Thus,  in  either  case  (a) 
when  Q{F)  =  0  or  (b)  when  Q{F)  >  0  it  follows  that  pi  =  P2  =  ps  =  0  on  [—1, 1].  g 


Remark.  Both  conditions  (a)  and  (b)  are  stated  in  terms  of  the  unknown  solution  F  of 
(Vc).  However,  since  the  shapes  of  the  functions  Fj  necessarily  fall  under  one  of  the  cases 
Cl  —  Cs,  it  is  possible  to  identify  hypotheses  on  the  boundary  values  F^  implying  that 
either  (a)  or  (b)  holds.  Specifically: 

(i)  Condition  (a)  is  equivalent  to  F~  =  F4,  j  =  1,2,3. 
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(ii)  If  Fj  <  F^,  Fj  >  F2  and  Fj  =  F^  then  Condition  (b)  is  satisfied. 

§3.4.  Existence  theorems  for  (Ve) 

We  conclude  this  Section  by  stating  two  existence  theorems  for  (Fe),  which  follow  from 
the  preceding  analysis.  The  first  applies  to  data  that  are  close. 

Theorem  3.7.  There  ia  a  r  >  Q  so  that  if  satisfy  (M)  and  \f'^  —  f~\  <  ihe 

boundary-value  problem  (F^)  has  a  solution.  The  parameter  r  may  depend  on  £. 

Proof.  Take  F(^)  =  f~  and  use  Lemma  3.6  part  (a)  and  Theorem  3.5.  | 

In  this  context  the  region  of  solvability  may  depend  on  e  and  thus  disintegrate  as  £  |  0. 
We  next  present  a  class  of  data  for  which  this  possibility  is  ruled  out. 

Theorem  3.8.  Let  /*  satisfy  (M)  and  /j“  </[*">  >  ft  ^  /T  —ft  •  £>  0  the 

boundary-value  problem  (F«)  has  a  solution. 

Proof.  A  continuation  eirgument  is  used.  With  the  parameter  /x  €  [0, 1]  let 

f-M  =  r ,  =  (/r + -  /D .  • 

Observe  that  =  0,  f'^{0)  =  f~,  /'’’(I)  =  /'*’  and 

frw  <  ffM .  .  /3'W  =  .  0  <  Z'  ^  1  ■ 


Define 


C  =  {n£  [0, 1]  :  (Fe)  with  data  f~ifi),  has  a  solution  }  . 
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Clezirly  0  €  C.  Furthermore: 

(i)  C  is  open. 

0  is  an  interior  point  of  C  (in  the  relative  topology),  by  Theorem  3.7.  If  now  €  (0,  l)nC 
let  Fft  be  the  corresponding  solution  of  (Ve)  and  consider  (3. 6-3. 7)  for  F  =  F^i.  By  virtue 
of  Lemma  3.6  part  (b)  and  Theorem  3.5  this  boundary-value  problem  has  a  solution. 
Therefore,  is  an  interior  point  of  C. 

(ii)  C  is  closed. 

This  is  an  immediate  consequence  of  Lemma  1.5. 

We  conclude  that  C  =  [0, 1]  and  in  particular  the  boundary- value  problem  {Vg)  with  data 
has  a  solution,  g 

Wc  emphasize  that  while  the  existence  of  solutions  part  requires  the  extra  assumptions 
/i"  <  /f  >  ft  1  fr  =  ft'i  the  fluid-dynamic  limit  part  (Theorem  2.1)  imposes  no 
restrictions  on  the  data,  aside  from  the  natural  requirement  that  /■*■  and  f~  are  strictly 
positive  Maxwellians. 

It  is  also  interesting  to  see  that  the  specied  data  in  Theorem  3.8  cannot  be  associated 
with  a  nontrivial  travelling  wave  solution  as  described  in  §0.  The  reeison  is  simple.  From 
(0.4)  we  see  a  travelling  wave  would  have  to  satisfy  the  limiting  conditions 

— •s(/i  +  /z  +  4/3)  +  {fi  —  f2)\  =  0  , 

\e=-oo 

-•s(/i  - /z)  +  (/i  + /z)  =0. 

0=—oo 

For  the  data  of  Theorem  3.8  this  would  imply  /j"  =  ft ■,  f2  ~  fti  fz  —  ft  which  yields 
an  everywhere  constant  solution  to  (0.4). 
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§4.  The  Carletnan  model 

The  Caxleman  model  of  the  kinetic  theory  of  gases  consists  of  the  system  of  semilinear 
hyperbolic  equations 

^f2  ~h  —  f2  f2\ 
dt  dx~ 

It  describes  a  system  of  two  kinds  of  particles  moving  with  velocities  +1  and  —1  respectively 
and  colliding  according  to  a  nonphysical  collision  rule.  As  a  model  in  kinetic  theory  it  has 
serious  defficiencies;  it  conserves  mass  but  not  momentum.  Nevertheless,  the  Carleman 
model  has  generated  some  interest,  because  it  enjoys  part  of  the  structure  of  the  Broadwell 
system  and  it  allows  for  rigorous  passage  to  the  hydrodynamic  limit  Kurtz  [11].  (See 
Platkowski  and  Illner  [12]  for  additional  references  and  a  discussion  of  relevant  results). 
In  this  section  we  explore  the  ideas  of  self-similar  hydrodynamic  limits  for  the  case  of  the 
Carleman  model. 

As  before,  a  modified  version  of  the  Carleman  is  considered,  admitting  solutions  of  the 
form  /(O  =  (/i(0>/2(0)>  with  ^  =  f ,  for  Riemann  data  /"  =  (/f./j"),  =  U\  J})- 

In  the  variable  ^  the  system  under  consideration  reads 


-{(.  - 1)/:  =  -QU) 


-(« + i)/i  =  - 


(4.2) 


(4.2) 


/(±1)  =  /*  ■ 

where  Q{f)  =  f 2  ~  fi  ■  The  data  are  assumed  positive  Maxwellians 


.  QU-)  =  Q(n  =  0- 


(4.4) 
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We  refer  to  this  problem  as  (T^e)- 


Our  goal  is  to  present  a  case  where  the  analysis  is  simple  and  can  be  carried  out  com¬ 
pletely,  with  no  further  restrictions  on  the  data.  Most  of  it  parallels  the  analysis  presented 
for  the  Broadwell  model,  and  our  exposition  will  be  sketchy  only  emphasizing  the  differ¬ 
ences. 


Definition.  The  pair  /(^)  =  hiO)  defined  on  [—1,1]  is  a  solution  of  (T^i),  if 

/i  €  C[-l,l]nC‘[-l,l),  /2  €  C[-l,l]nC‘(-l,l]  satisfy  for  ^1,  6  €[-1,1]  the  integral 
equations 


-(^+l)/2 


(4.5) 


and  the  boundary  conditions  /j(±l)  =  /j^,  j  =  1,2. 


Theorem  4.1.  For  each  pair  of  positive  Maxwellian  data  /“,  /■*■  and  for  each  e  >  0 
the  boundary  value  problem  (TZ,)  has  a  solution.  The  solution  f  is  Holder  continuous  on 
[— 1, 1]  with  exponent  Og  {Lipshitz  continuous  for  e  <  eo(/^))- 


Let  now  /*  be  fixed  and  consider  a  family  {f‘}e>o  of  solutions  to  {TZg).  The  functions 
/*  are  extended  to  (—00,00)  by  setting  /*  =  f~  on  (— 00,— 1)  and  /'  =  on  (l,oo). 
The  extended  functions  are  called  again  /'  and  satisfy  (4.5)  for  ^i,  ^2  €  (—00,00). 


Theorem  4.2.  Let  {/'}e>o  be  a  family  of  extended  solutions  of  (TZg)  taking  fixed,  positive 
Maxwellian  data  /^.  There  exists  a  subsequence  {/'"  }  with  ♦  0  and  a  positive,  bounded 
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function  p{^)  such  that 


(/!"(«)  ■*  5  (m«)  ^ 

The  function  p(^)  is  a  weak  solution  of 

'2f-  X  <  0 

dtp  =  Q  ,  p(x,0)  =  < 

.2/+  x>0 

that  is 

(2f-  X  <  0,  t  >  0 

P(j)=< 

I  2/+  X  >  0,  <  >  0 

Proof.  The  strategy  for  proving  Theorems  4.1  and  4.2  is  analogous  to  the  Broadwell  case. 
The  proof  is  decomposed  in  three  parts. 


Part  1.  A-priori  estimates 

Let  /  be  a  solution  of  (T^e).  Using  (4.5)  it  follows  that  Q{f{-1))  =  Q(/(+l))  =  0.  Since 
Q{f)  satisfies  the  differential  equation 


(4.7) 


either  Q{f)  =  0  on  (—1,1)  or  it  never  vanishes.  Hence,  the  shapes  of  the  components 
/i ,  /2  are  restricted  as  follows: 

^1-  Qif)  >  0  on  (—1, 1).  Then  /i  and  /2  are  both  increasing  on  (—1, 1). 

Qt/)  <  0  (“L  !)•  Then  f\  and  /j  are  both  decreasing  on  (—1, 1). 

^3:  Q{f)  =  0  on  (—1,1).  Then  /i  and  /j  are  constant  functions. 
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As  a  consequence,  for  ^  €  [~1,  Ij 


(4.8) 


niin{/j  ,  /+)  < /j  <  max{/j  , /+ }  ,  j  =  1,2 

TVfi  =  \ff-f;-\  ,  i  =  l,2. 

Next,  Eui  argument  similar  to  Lemma  1.3  shows  that  there  are  expenents  pe,  1  <  Pc  <  oo, 
and  positive  constants  Cg  such  that 


II/:- lip.  <  c-c 


i  =  i,2. 


(4.9) 


The  exponent  pt  depends  on  /*  and  e,  it  is  increasing  as  £  decreases,  and  there  is  So  = 
£o(f^)  >  0  so  that  for  £  <  £o  the  exponent  pe  =  oo.  It  follows  from  (4.9)  that  solutions  of 
(T^e)  are  Holder  continuous  (Lipschitz  continuous  for  £  <  £o). 


Part  2.  The  fluid  dynamic  limit 

Let  {/*}e>o  be  a  family  of  extended  solutions  to  (T^g)  corresponding  to  fixed,  positive 
Maxwellian  data  /^.  By  (4.8)  and  Helly’s  theorem,  there  exists  a  subsequence  {/'"}  with 
£„  — >  0  and  a  function  /  =  (/i,/2)  with  positive,  bounded  components  /i,  /2  of  bounded 
variation  such  that  /'"  — >  /  pointwise  on  (—00,00).  Clearly,  /(^)  =  f~  for  ^  6  (—00,  — 1] 
and  /(O  =  for  (  €  [l,oo). 

The  members  of  the  sequence  satisfy  the  identities 


(/f"  +  /!")'  +  {/;•  -fi")'  =  0 


+  f!"  Infi' ) ’  +  (/f"  = 

=  -  A'"')  • 


(4.10) 


(4.11) 
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Identity  (4.11)  is  obtained  upon  multiplying  the  first  equation  (4.2)  by  (1  +  infl”),  the 
second  by  (1  +  inf 2')  and  adding  the  resulting  equations. 

Fix  now  ^1,  ^26  [—1, 1]  and  combine  (4.11)  with  the  uniform  bounds  (4.8)  to  obtain 

-EnC  <  -  f'  (/f"’  -  /I"") 

<-i  r  (f;~ 

and  upon  passing  to  the  limit  e„  j,  0 


=  0. 


We  conclude  fi  =  for  a.e.  ^  €  (-00,00). 

Next,  set  p  =  f\  f 2-  Passing  to  the  limit  in  (4.10)  we  see  that  —^p'  =  0  weakly. 

Therefore  p{~)  is  a  weak  solujtion  of 


5«/9  =  0  ,  p(x,0)  = 


2/  X  <  0 


2/+  X  >  0 


(4.12) 


The  solution  of  (4.12)  is  given  by 


(  2f~  X  <  0,  t  >  0  , 


P  =  p(-^)={ 


[2f^  X  >  0,  0  0  . 


Remark.  It  follows  from  (4.6)  suid  (4.8)  that  for  any  1  <  p  <  00 


As  a  matter  of  fact,  since  the  limiting  function  p  is  characterized  as  the  unique  solution  of 
(4.12),  the  above  convergence  holds  over  the  entire  family  {/'}e>o  as  e  0+. 
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Part  3.  Existence  theory  for  (T^e). 

First,  we  introduce  some  notation.  Let  h{f)  =  (— 2/i,2/2)  and  define  the  matrices 


'-(«-!)  0  ■ 

■-2/1 

2/2  ■ 

0  -(«  +  !). 

»  ^u)  — 

.  2/1 

-2/2. 

Then  Q{f)  =  ■  f  and  (4.2)  may  be  written  as 


=  ^BU)f  ■ 


(4.13) 


Let  F  =  {F\ ,  F2)  be  a  given  solution  of  (T^e)  defined  on  [—1, 1]  and  taking  boundary  val¬ 
ues  F(±l)  =  F^  positive,  Maxwellian  states.  By  (4.9)  the  solution  F  is, Holder  continuous. 
Any  nearby  solution  /  may  be  written  as  f  =  F  +  p,  with  the  perturbation  p  =  (pi,P2)  a 
continuous  function  satisfying 


A(Op'  =  jB(F{0)p  +  ^B(p)p 

p(±l)  =  p* 

and  p*  =  /^  —  F*.  Since  /*  and  F^  are  Maxwellians 


(4.14) 


b(F^)-p^  +  Q{p^)  =  0.  (4.15) 

The  first  step  in  proving  Theorem  4.1  is  to  show  the  following  analog  of  Theorem  3.5. 
Theorem  4.3.  If  the  only  continuous  solution  of  the  homogeneous  boundary-value  problem 


/!({)?'  =lB(f(«))p , 


(4.16) 


p(±l)=0, 

is  the  trivial  solution  p  =  0,  then  there  exists  a  positive  constant  r  such  that  for  any  data 
satisfying  (4.15)  and  |p"^|  4-  |p~|  <  r  the  nonlinear  boundary-value  (4.14)  has  a  solution. 
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Proof.  The  proof  is  long,  but  it  follows  the  same  steps  as  the  proof  of  Theorem  3.5  being 
just  easier.  We  only  give  a  rough  sketch  here. 

First,  one  analyzes  the  singular  linear  boundary-value  problem 


X(e)p'  =  iB(f«))p  +  lG(0, 

e  £ 

K±i)  =  p* , 


(4.17) 


with  subject  to  the  restrictions 


6(F+)  •  -b  pi(-bl)  =  -b{F-)  •  p-  +  P2(-1)  =  0  .  (4.18) 


It  turns  out  that  (4.17)  admits  a  unique  solution  for  any  data  p*  subject  to  (4.18)  if  and 
only  if  (4.16)  possesses  only  the  trivial  solution  p  =  0.  Moreover,  in  that  case  the  solution 
p  of  (4.17)  can  be  estimated  by 

l|p||<C(|p+l  +  |p-|  +  ||Gi|),  (4.19) 


where  ||  •  ||  stands  for  the  sup-norm  on  [—1, 1]. 

Next,  one  applies  the  contraction  mapping  principle  to  a  map  carrying  P  to  the  solution 
of  (4.17)  with  G  =  ^B(P)P.  Given  the  estimates  (4.19),  the  proof  of  this  part  is  identical 
to  the  proof  of  Theorem  3.5.  The  resulting  fixed  point  is  the  solution  p  of  (4.14).  | 

The  second  step  is  to  show: 

Lemma  4.4.  The  only  continuous  solution  of  the  homogeneous  boundary- value  problem 
(4.16)  is  the  trivial  p  =  0. 
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Proof.  Suppose  that  p  =  (pi,p2)  is  a  nontrivial  eigensolution  of 


[(1  -  Opi  -  (^  +  1)P2]'  +  (pi  +  P2)  =  0  , 
and  integration  over  any  [a,  6]  C  [-1, 1]  yields 

(1  -  6)pi(6)  -  (1  -  a)pi(a)  -  (6  +  1)p2(6)  +  (a  +  l)p2(a)  +  /  (pi+p2)d^  =  0.  (4.21) 

J  a 

Consider  the  curve  7  :  P  =  (pi(0,P2(0),  ^  €  [-1, 1],  that  describes  the  solution  in  the 
state  plane  pi  -  p2.  Since  (0,0)  is  an  equilibrium,  7  only  meets  the  origin  at  the  singular 
points  ^  =  ±1.  The  vector  field  in  (4.20)  vanishes  only  on  the  lines  Fi(4)pi  -  Fzi^pt  =  0 
lying  in  the  first  and  third  quadrants.  This,  together  with  the  form  of  (4.20),  implies  that  7 
may  only  cross  the  positive  p2-axis  going  from  the  second  to  the  first  quadrant,  the  positive 
Pi -axis  going  from  the  first  to  the  third,  the  negative  p2-axis  from  the  fourth  to  the  third, 
and  finally  the  negative  pi-axis  from  the  third  to  the  second.  There  are  the  following,  not 
mutually  exclusive,  possibilities  to  consider:  (i)  7  crosses  the  positive  p2-ajcis,  (ii)  7  crosses 
the  negative  p2-axis,  and  (iii)  7  starts  and  concludes  at  the  origin  without  crossing  the 
P2-axis. 

Case  (i):  Suppose  that  p  crosses  the  positive  p2 -axis  at  ^  =  a.  Thenpi(a)  =  0,  P2(a)  >  0. 
Also,  the  curve  7  will  either  cross  the  positive  pi-axis  at  some  ^  =  b,  or  it  will  wander 
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around  in  the  first  quadrant  concluding  on  the  origin  at  ^  =  6  =  +1.  In  either  case 
pi{b)  >  0,  Piib)  =  0  and  pi(0  >  0)  P2{0  >  0  a  <  ^  <  b.  Identity  (4.21)  then  leads  to 
a  contradiction  and  that  excludes  case  (i).  Case  (ii)  is  similarly  excluded  by  applying  the 
same  argument  to  (— p),  which  is  again  a  solution  of  (4.20). 

We  conclude  that  7  starts  and  concludes  at  the  origin  without  crossing  the  pi-axis. 
Without  loss  of  generality  we  may  assume  that  the  curve  7  lies  on  the  right  half  plane; 
otherwise  replace  p  by  (— p)  as  before.  Since  pi(  —  l)  =  pi(+l)  =  0,  the  mean  value  theorem 
and  (4.20)  imply  there  exists  6  €  (  —  1, 1)  such  that  Fi{d)p\{6)  =  F2{9)p2{9).  At  9  the  curve 
7  lies  in  the  first  quadrant.  In  view  of  the  form  of  the  vector  field  in  (4.20),  the  curve  7 
starts  at  the  origin  at  ^  =  —1,  lies  in  the  fix'st  quadrant  thereafter  up  to  at  least  ^  =  9, 
and  then  either  escapes  to  the  second  quadrant  at  some  point  ^  =  6,  or  it  concludes  at 
the  origin  at  ^  i  =  +1.  In  either  case,  pi(6)  >  0,  P2(b)  =  0  and  pi((f)  >  0,  P2(0  >  ^ 
for  -1  <  ^  <  6.  Applying  (4.21)  with  a  =  —1  and  b  as  above,  we  arrive  again  at  a 
contradiction.  Therefore,  p  =  0  on  (—1, 1].  | 

The  last  step  in  proving  Theorem  4.1  is  to  use  a  continuation  argument.  Let  f~  = 
(/r>/2~)i  f'^  =  positive  Meixwellians,  that  is  7/2” i/i*” 7/2^  >  0  /f  = 

=  ft  ■  For  p  €  [0, 1],  set 

r  (p)  =  r  7  r  (p)  =  r + p(r  -n, 

which  are  also  positive  Maxwellians.  Consider  the  set 

C  =  {p  €  [0, 1]  :  (T^e)  has  a  solution  taking  data  /“(p),  ■ 
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Then  0  €  C  and,  by  virtue  of  Theorem  4.3  and  Lemma  4.4,  C  is  open. 

We  proceed  to  show  that  C  is  closed.  Let  {^„}  be  a  sequence  of  points  in  C  such 
that  /in  — ♦  /io-  Suppose  that  are  the  corresponding  solutions  of  (T^i)  admitting  data 
/■•■(/in)  and  /“(/in)-  On  account  of  (4.S)  and  (4.9),  the  sequence  {/''"}  is  precompact 
in  some  Holder  space.  Therefore  a  subsequence  converges  to  a  limiting  function  ,  and 
a  straightforweird  argument  shows  that  /^®  is  a  solution  of  (J^t)  taking  data  /''"(/io)  and 
/"(/<o)-  Thus  C  is  closed. 

We  conclude  that  C  =  [0, 1]  and  the  proof  is  complete,  g 

§5.  Appendix:  A  linear  singular  equation. 

We  record  here  certain  properties  of  the  equation 

-«p'  +  <iK)p  =  ft({)  (5-1) 

that  are  used  in  the  text  to  study  the  behavior  of  solutions  near  singular  points.  The 
equation  is  taken  in  the  interval  0  <  <f  <  a  with  a  <  1.  The  objective  is  to  investigate 
regularity  properties  for  solutions,  and  to  study  the  map  that  carries  h  to  p. 

Throughout,  the  functions  (f>  and  h  are  cissumed  at  least  continuous  on  [0,  a]  with  0 
subject  to  the  restrictions; 

„0)>0,  r|m^i.c<co.  (H.) 

Jo  s 

The  latter  assumption  is,  for  instance,  satisfied  if  the  function  <j)  is  Holder  continuous. 
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Use  of  an  integrating  factor  shows  the  solution  of  (5.1)  on  (0,  a]  is  given  by 

r  his) 

P(0  =  P(°)  /  - -  eJ‘  <  ^  ds . 

Next,  certain  convenient  expressions  for  the  solution  are  derived.  Observe  that 

J(()  :=  e“^  ^ 


where  u)  is  defined  in  terms  of  the  convergent  (due  to  (Hi))  integral 


-i: 


m  -  m 


On  account  of  (5.3)  and  (5.4),  p  may  be  written  in  the  form 

p(^)  /  h{s)  ds  . 

a 

An  alternative  formula  for  p  is  obtained  by  combining  (5.5)  with  the  identity 
/  ds  = 


_1_  a-<^(0)  g-w(a)  _  ^-0(0)  r  ^-,^(0)  g-u>(s)  <A(^)  -  <^(0) 

<^(0) 


(5.4) 


It  reads 


^  ^<^(0)  f  ^(■^)  ^(0)  „-<»(0)  ^u>(C)-w(5) 

Jf  5 


M21  c'-AtO)  *«^(0-U/(5) 


u  s 


5  '  e 


Formulcis  (5.5)  and  (5.6)  serve  as  starting  points  to  produce  various  estimates  for  solutions 
of  (5.1). 
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Lemma  5.1.  Let  <f>  and  h  be  continuous  functions  on  [0,  a]  satisfying  (i?i).  Then 


limp(0  = 

p  is  continuous  on  [0,  a]  and  satisfies  the  bound 

+  sup  |ft(s)e-“<-)|. 

0<s<a 


(5.7) 


(5.8) 


Proof.  Hypothesis  (Hi)  implies  that  u;(^) 
h(s)  ds 


lim 

€-0 


^-Mo) 


4  0  as  ^  0.  Then  L’Hopital’s  rule  shows 

=M0)M0), 


which  together  with  (5.5)  implies  (5,7). 

Ne.xt,  combining  the  estimate 

|p(0e-‘*’(^M  <  |p(«)e"‘"^“M+  sup  |/»(s)e-‘*'<^>|(e^^°)  f  , 

0<3<a 

obtained  directly  from  (5.5),  with  the  formula 


^.(0)  /%-*(»)-.  ds  =  ^(1  -«/„)*"») 

we  deduce  (5.8).  | 


The  next  task  is  to  relate  the  modulus  of  continuity  of  <f)  and  /i  at  ^  =  0  with  the  modulus 
of  continuity  of  the  solution.  To  this  end,  eissume  that  for  some  0  <  /?  <  1  and  0  <  7  <  1 
the  functions  ^  and  h  satisfy 


|^(O-fl^(0)l<<<^> 

|/i(0 - M0)l  <<h>  c 


(H2) 
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on  [0,  a];  here,  the  symbols  <  <f>  >  and  <  h  >  stand  for  the  Holder  constants.  We  show: 


Lemma  5.2.  Let  (f>  and  h  satisfy  (Hi  —  H2).  If  =  7  =  1  and  (^(0)  >  1  then  p  is  Lipshitz 
and  satisfies  the  bound 

|p(0  -  P(C)|  <  C  {  »-«<»)  |p(a)  -  ^  I  +  <  ft  >  +  I  ^  I )  If  -  Cl  e ,  C  .  [0,  a] .  (6.9) 

Otherwise,  p  is  Holder  continuous  with  exponent  any  0  <  a  <  oq  :=  min{/?,7,  ^(0)}  and 
satisfies  for  ^ ,  C  f  [0i  a] 

|P(0-P(C)|  <  C  -  -^1+  <  ^  >  l^-Cr-  (5.10) 

The  constant  C  depends  on  <  <^  >,  /?,  7,  a  and  ^^(0),  but  is  independent  of  a  in  the  region 

0  <  a  <  1. 

Proof.  In  what  follows  C  will  stand  for  a  generic  contant  exhibiting  the  aforementioned 
dependence.  First,  the  terms  in  the  right  hand  side  of  (5.6)  are  estimated.  In  view  of  (5.4) 
and  (H2) 

Kf)  - ‘.'Wl  <  <  •^  >  I =  {5-11) 

Since  a  <  1,  (5.6)  and  (5.11)  imply  for  0  <  ^  <  a 


\piO  - 


m 

m 


<|p(„)_M|e^  (!)'«”> +  <A>  (f*<")  <i^) 

</>(0)'  a 

1^1 
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where 


(5.12) 


Set  now  Qo  =  niin{/?,7,(^(0)}  >  0  and  note  that  0  <  ao  <  1.  Using  the  formula  for  the 
evaluation  of  the  integrals 


^0(0)  J'  ^a-l 


(  \ln^\ 


if  a  =  <^(0) 


«  1]^ 


(5.13) 


we  conclude 


h{0) 


ifao  =  <^(0) 


if  oo  <  ^(0) 


for  ^  e  [0,  a]. 


(5.14) 


The  control  on  the  modulus  of  continuity  of  p(^)  at  ^  =  0  provides  bounds  on  the 
derivative.  To  show  that,  rewrite  (5.1)  as 


(p'io  =  mpio  -  HO 

=  ^  (<>(«)  -  <«0))  +  ^(0)  (p(0  -  -  (HO  -  MO)) 

+  (<»(0  -  <>(0))  (p(0  -  (5.15) 


and  use  {H2)  and  (5.14)  to  obtain 


IKI 


if  Oo  =  <^(0) 
ifQO<<^(0)  ’ 


for  ^  €  (0,  a  . 


(5.16) 


We  emphasize  that  the  constant  C  in  (5.16)  depends  solely  upon  <  <f>  >,  /3,  -y  and  upper 
bounds  for  <^(0).  All  other  dependences  are  incorporated  in  the  constant  A  in  (5.12). 

We  can  now  complete  the  proof  of  the  lemma: 
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(i)  If  yS  =  7  =  1  <  ^(0)  then  ao  =  1  <  ^(0)  and  (5.9)  follows  from  (5.16). 

(ii)  In  all  other  cases,  fix  any  0  <  a  <  q;o  <  1  and  let  (e  (0,  a].  Using  Holder’s  inequality, 
we  deduce  from  (5.16) 


M()-p(0\  =  \  J%'U)ds\ 

<1  f\p'U)\'^ds\'-“]i-cr 

i-an 

f  s  (1  +  \lns\)'^^  ds  1^  —  Cr  •  (5-17) 

0 

Since  <  1,  the  last  integral  is  finite  and  (5.10)  follows  from  (5.17)  and  (5.7).  | 
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